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A novel type of symmetry breaking is reported for the two-channel Kondo lattice where conduction 
electrons have spin and orbital (channel) degrees of freedom. Using the continuous-time quantum 
Monte Carlo and the dynamical mean-field theory, a spontaneous breaking of the orbital symmetry 
is observed. The tiny breakdown of orbital occupation number, however, vanishes if the conduction 
electrons have the particle-hole symmetry. The proper order parameter instead is identified as a 
composite quantity representing the orbital-selective Kondo effect. The single-particle spectrum 
of the selected orbital shows insulating property, while the other orbital behaves as Fermi liquid. 
This composite order is the first example of odd-frequency order other than off-diagonal order 
(superconductivity), and is a candidate of hidden order in /-electron systems. 



Localized / electrons with even number per site can 
have a non-Kramers doublet ground state. If these 
/ electrons couple with conduction electrons with or- 
bital degeneracy, a two-channel Kondo system can be 
realized Recently two-channel Kondo systems have 
been attracting renewed attention, partly due to discov- 
ery of new Fa doublet compounds such as PrIr2Zn2o 
and FrVTAlon 0-13] ■ Mysterious behaviors observed in 
PrAg2ln and URu2Si2 and their impurity coun- 
terparts such as Pr2.Lai_2,Ag2ln [3] and Ua;Thi_a;Ru2Si2 
have remained long-standing puzzles, and the two- 
channel Kondo effect has been suspected as an origin. 

The two-channel Kondo iinpurity has the residual en- 
tropy in the ground state Then the two-channel 
Kondo lattice (2ch KL) should undergo some phase tran- 
sition to reach the ground state from high-temperature 
phase. Possible orderings in the 2ch KL have been stud- 
ied by the dynamical mean- field theory (DMFT) [13- 
[11], and by other methods [l^l- Among possible ground 
states, the simplest is a magnetically ordered state [lo| . 
Superconductivity with odd frequency pairing has also 
been discussed for the 2ch KL 11] and a related model 
called the SU(2)xSU(2) Anderson lattice [3. On the 
other hand, an exact diagonalization study, using a small 
discretized reservoir in the DMFT, has found a channel 



symmetry breaking 12] • In the effective impurity prob- 
lem of the DMFT, a minute energy scale appears as- 
sociated with the Kondo effect. Then highly accurate 
treatment of the reservoir is required for numerical ap- 
proach. Recently, a Monte-Carlo scheme using the con- 
tinuous imaginary time (CT-QMC) [isl] has proven pow- 
erful enough in sol ving the effective impurity problem at 
finite temperature [id [l7j. 

In this letter, by combining DMFT and CT-QMC, we 
report on a new kind of electronic order in the 2ch KL at 
and near half-filling of two equivalent conduction bands. 
Although the symmetry of channel occupation is broken, 
the corresponding order parameter is tiny and vanishes 
in the case with the particle-hole symmetry. The proper 
order parameter is a composite quantity that becomes 
order of unity. To our knowledge, this is the first example 



of the composite order parameter other than off-diagonal 
order (superconductivity) [13, 3-23]. It is shown that 
the present composite order can be interpreted also as an 
odd-frequency order. 

The simplest description of periodic non-Kramers dou- 
blets coupled with conduction electrons is given by the 
2ch KL [iMl: 



n 



E 



(1) 



where Ckatr is annihilation operator of the conduction 
electron with momentum fc, channel a = 1, 2 and pseudo- 
spin (7 =ti i- The localized pseudo-spin at site i is de- 
scribed by the spin 1/2 operator Si, while Sda represents 
the conduction pseudo-spin of channel a at site i. For the 
moment, we call the pseudo-spin simply as 'spin'. Near 
the end of the letter, we return to physical consequence of 
the non-Kramers doublet described by the pseudo-spin. 

In our DMFT, the semi-circular density of states 
is taken for conduction electrons defined by Pq{s) — 
(2/7r)-y/l — {e/DY with the nesting property Sk+Q — 
—Ek where Q is a staggered ordering vector. We take 
D = 1 as a unit of energy. At half filling of conduction 
bands, the system has a particle-hole symmetry. We have 
employed the Pade approximation for analytic continua- 
tion to real frequencies. 

The phase diagram of the 2ch KL near = 2 (half fill- 
ing) is shown in FIG. [T] where ric is a number of conduc- 
tion electrons per site. Figure [Ija) shows the transition 
temperature of antiferro(AF)-spin order that has been 
discussed in Ref. [ll[. With lower transition tempera- 
ture at half-filling, a ferro(F)-channel order also emerges 
as shown in FIGllIb), as has been pointed out in Ref. 
[13] ■ In the present case, the F-channel phase is masked 
by the AF-spin phase around half-filling because of the 
lower transition temperature. In real systems, however, 
one may encounter possible suppression of the AF-spin 
order caused, e.g., by geometrical frustration, or substan- 
tial next-nearest neighbor hopping. In such a case, the 
F-channel order may be stabilized. Because of its unique 
and interesting property to be described below, this let- 
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FIG. 1: (color online) Phase diagram of the 2ch KL near half 
filling for (a) AF-spin and (b) F-channel ordered phases. In 
(b), we neglect the AF-spin order for all ric. The lower panels 
show spin and channel moments close to the ground state, as 
defined by eqs. and ((S}. The full moment is normalized 
to unity. 
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FIG. 2: (color online) (a) Temperature dependence of lo- 
cal correlation functions, and (b) schematic picture of the 
F-channel phase. The arrows on the thin lines show conduc- 
tion electrons, and the shaded ovals show the Kondo singlets 
centered on each lattice site, (c) Single-particle spectra of con- 
duction electrons with a = 1 and 2 in this phase are shown 
in (cl) and (c2), respectively. 



ter concentrates mainly on the F-channel phase from now 
on. 

Let us define the spin and channel moments: 



a 



(2) 
(3) 



where is the local number operator of conduction 
electrons with channel a and spin a. The lower panels of 
FIG. [1] show calculated results of the moments. The spin 
moment ai ric — 2 in AF-spin phase becomes maximum 
with highest transition temperature Tgpin, and gradually 
decreases away from half filling. In the case of F-channel 
phase, on the contrary, the channel moment mchan be- 
comes finite only away from half filling, and remains tiny. 
Even though the transition temperature Tchan takes the 
maximum at ric — 2, we observe mchan = 0. Hence, the 
channel moment is not a proper order parameter. 

Let us identify the proper order parameter in the F- 
channel phase. It has been found that the double occu- 
pancy {riafriai) in the ordered phase becomes different 



between a = 1 and a = 2 [l^]. We propose, however 



that the fundamental order parameter involves the local- 
ized spin, and hence the Kondo effect. The small differ- 
ence in double occupation comes out as a consequence of 
the fundamental order parameter. We shall demonstrate 
that the proper order parameter leads to identifying an 
odd-frequency order. As shown in FIG. ^a), the local 
spin correlations {Si-Sda) becomes different below Tchan- 
Namely each localized spin tends to form the Kondo sin- 
glet selectively with one of the two conduction bands. 
The order parameter is hence given by 



^ = {Si - (Sctl - Sci2)), 



(4) 



which is independent of site index i. Since 5* grows con- 
tinuously below the transition temperature, the phase 
transition is of second order. Note that the order param- 
eter is a composite quantity, and cannot be described 
by a one-body mean field such as (5'^) or {sdi — Sci2)- 

Real-space image of the electronic state is illustrated 
in FIG. [ll^b). For channel a = 1, the effective Kondo 
coupling tends to zero, while for a = 2 the coupling tends 
to infinity. Thus the F-channel phase is the mixture of 
weak- and strong-coupling limits depending on channels. 
This state therefore cannot be accessible by perturbation 
theory from either limit. 

The peculiar character of the F-channel phase appears 
also in the single-particle spectrum. We have derived the 
single-particle spectrum explicitly from the imaginary 
part of the Green function. Since the self-energy is local 
in the DMFT, the wave- vector enters only through e^- 
We introduce the parameter k defined by Ek = cos k, 
and visualize the spectrum as if the system were in one 
dimension. Accordingly the single-particle spectrum is 
written as A{k,u!). Figures[5I^cl) and (c2) show the spec- 
tra of conduction electrons with a = 1 and a = 2, respec- 
tively. The spectrum of the channel a — 1 displays the 
Fermi-liquid behavior. Here the mass enhancement fac- 
tor is estimated as m*/m ~ 1.95 from analysis of the 
self-energy. As shown in FIG.[2Jc2), on the contrary, an- 
other channel a = 2 acquires the insulating character. 
The spectrum is almost the same as that of the ordinary 
Kondo insulator. Thus, the F-channel phase consists of a 
Fermi liquid with a = 1 plus Kondo insulator with a = 2. 
Hence the phase transition at T = Tchan is regarded as 
an orbital-selective metal-insulator transition. The rea- 
son for metal and insulating characters of each channel 
is understood easily from FIG.dJb). 
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FIG. 3: (color online) (a) Spin- and (b) channel-dependent 
differences of local Green functions defined in eqs. ([Sj and 
The differences of complex self energies defined in eqs. 
(dJ and ((9)1 are shown in (c) and (d). 



We proceed to show that the F-channel ordered state 
can be interpreted as an odd-frequency order, which is a 
consequence of the composite nature of \E' given by (H)) . 
Let us take the case ric = 2 with particle-hole symmetry, 
and compare the F-channel order to the AF-spin order. 
We take the difference of local Green functions and in- 
troduce the quantities: 

Aapm(r)=^[G^^(r)-G^,(r)], (5) 

a 

AGcha„(T) = ^[Gi,(t) - G2.(t)], (6) 
a 

where G^„(t) = -(r^Q„<x(r)ct^^). In the AF-spin 
phase, rngpin = AGspin(O) becomes nonzero, and serves 
as the order parameter. In the F-channel phase, on the 
other hand, we have seen in FIGHJb) that the apparent 
order parameter rrichan = AGchan (0) is tiny and vanishes 
at ric = 2. 

Figure |3] (a) and (b) show the r-dependence of these 
quantities for the case of tt-c = 2. It can be seen that 
AGspin (t) in the AF-spin phase is an even function of 
T, while AGchan(''') in the F-channel phase is an odd 
function. Hence the derivative dAGchan/dr at r = 0, 
rather than mchan, is the proper quantity to describe the 
F-channel order. The explicit form is given by 



aAG, 



cha 



dr 



J 



where N is the number of sites. The second term in the 
right-hand side (RHS) corresponds precisely to the order 



parameter ^ introduced by ([4]). On the other hand, the 
kinetic energy given by the first term also becomes depen- 
dent on channel. However, this quantity is not suitable 
as a proper order parameter, since it does not involve 
localized spins. 

The low-energy limit of the self-energy corresponds to 
a one-body mean field at the Fermi level. Then we intro- 
duce the difference of self-energies as 

AE3pin(ie„) = EP"t(i£") " ^U^^n)], (8) 

a 

AEchan (i£«) = ^[i;io-(ie„) - i;2<T(i£«)], (9) 

(7 

where e„ = (2n + l)7rT is the fermionic Matsubara fre- 
quency. Figure |3] (c) and (d) show the computed re- 
sults for AI](ie„). For AF-spin phase, AEgpin is real and 
even function of Matsubara frequency e„. For F-channel 
phase, on the other hand, AEchan is pure imaginary and 
odd function of e„. Hence, the F-channel order is identi- 
fied as an odd-frequency order also from the self-energy. 
It should be noted that the present clear-cut classification 
applies only to systems with the particle-hole symmetry. 
Otherwise, both orders involve finite components of the 
opposite parity, and are neither completely even nor odd. 
We emphasize that the concept of the composite order is 
valid even without the particle-hole symmetry. 

Let us now discuss the AF-spin and F-channel orders 
in terms of susceptibilities. We first introduce the two- 
particle Green function x^^^., (ie„,ie„/;iz^m) where = 
2rmTT is an bosonic Matsubara frequency. The AF-spin 
and F-channel orders correspond to wave vectors q = Q 
and q = 0, respectively. The relevant two- particle Green 
functions for our purpose are obtained by appropriate 
linear combinations with respect to channel and spin in- 
dices in Xqaa'i ^ud are written as x^Tn(i^n' i^n'! i^™) ^^'^ 
Xchan(iSn' i^n' ! i'^m)- Then the following even- and odd- 
frequency susceptibilities are obtained: 



Xcvcn — Q ^ ^ X(i£n; i^n^ ; Q) ; 
nn' 



(10) 



Xodd 



13 



where we have suppressed the indices such as AF and 
spin. The even part Xcvon coincides with the conven- 
tional susceptibility, and describes fiuctuations of spin 
and channel moments. The odd function g{en) in eq. 
([TT|) abstracts the odd-frequency part of the two-particle 
Green function 



11, 2 



In order to determine the appropriate g{en), we intro- 
duce an operator O by Q 



N ^ 



' dr 



dCi2a(T) 



W2cr 



dr 



(12) 
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FIG. 4: (color online) Temperature dependence of 1/xcvcn 
and 1/xodd for (a) AF-spin and (b) F-channel orders. Phase 
transitions are signaled by the zero crossing of inverse suscep- 
tibilities. 



The expectation value (O) coincides with the derivative 
([7]). With use of O, one can prove the following relation: 



N / (C'(T)C't)dT 



where convergence factors enter in frequency summation. 
The first term in the RHS must diverge at the phase tran- 
sition. It can be shown that replacing e„ by (?(£„) leads to 
the divergence at the same transition temperature, pro- 
vided g{en) has finite overlap with e„ as a vector in fre- 
quency space. For simpler calculation, we want to avoid 
the convergence factor, and choose g{en) = tanh(£„/£)) 
in eq. (fTTjl . Because of the minus sign in the first term 
of the RHS of ([T3|) . Xodd becomes negative infinity at the 
transition temperature. 

Figure m shows temperature dependence of inverse sus- 
ceptibilities in the particle-hole symmetric case — 2. 
For the AF-spin order as shown in FIG. |4] (a), Xcvcn 

di- 
verges at Tspin ~ 0.082 indicating the spin order below 
^spinj while Xodd docs not show any anomaly. On the 
other hand, for the F-channel order as shown in FIG. 
Sl(b), Xcvcn does not signal any instability. This corre- 
sponds to mchan = in the F-channel phase. However, 
Xodd shows drastic temperature dependence. With de- 
creasing temperature, 1/xodd becomes zero at Tchan — 
0.027 from negative side. This means that Xodd defined 
by pT|) indeed signals the divergent fluctuation of com- 
posite order parameter 5*. 

Finally, we discuss possible relevance of our results in 
understanding real /-electron systems. If we interpret 
the pseudo-spin as describing a non-Kramers doublet, 
the channel index a = 1 is translated into the real up 
spin, and a = 2 into the down spin. Here, the choice of 



up and down spins can of course be reversed. Then the 
AF spin phase corresponds to AF orbital order, while 
the breakdown of channel symmetry is translated into 
that of the real-spin symmetry. Namely, the time re- 
versal symmetry is spontaneously broken in F-channel 
phase of the non-Kramers doublet lattice. Without the 
particle-hole symmetry, the present order gives rise to a 
tiny magnetic moment corresponding to mchan in FIGH] 
(b). However, the proper order parameter is given by 
which actually represents an octupole moment involving 
conduction electrons. Namely, spatial extent of up and 
down spins become different from each other, since con- 
duction electrons with up spin forms a Fermi liquid, while 
those with down spin tend to form orbital Kondo singlets 
with localized non-Kramers doublets. The resultant spin 
imbalance sums up to zero, but gives rise to extended 
magnetic octupoles centered at each site. 

Since 4* is independent of sites, any diffraction mea- 
surement cannot detect the present order. Nevertheless, 
transport properties signal clear phase transition involv- 
ing metal-insulator transition of down spins. The ordered 
phase is a kind of half-metal since down spins are insu- 
lating. Moreover, the size of the specific heat anomaly 
should be substantial since is of the order of unity 
as shown in FIGlSl^a). From these characteristics, it is 
tempting to consider ^E" as a candidate of the hidden or- 
der parameter in URu2Si2. Of course, much more refine- 
ment of the model is necessary before we can push this 
suggestion further. 

In conclusion, we have demonstrated a novel type of 
symmetry breaking between equivalent channels in the 
2ch KL. The ordered state is characterized by a two- 
body correlation function, or composite order parameter. 
In a general case without the particle-hole symmetry in 
conduction bands, tiny channel moment also arises as a 
deceptive order parameter. The composite order is re- 
garded as an odd-frequency ordering, which is the first 
explicit example other than superconductivity. 
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